Fully self-consistent ion-drag force calculations for dust in collisional plasmas with 

external electric field 
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The ion-drag force on a spherical dust particle immersed in a flowing plasma with external electric 
field is self-consistently calculated using the Particle In Cell code SCEPTIC in the entire range of 
charge-exchange coUisionality. Our results, not based on questionable approximations, extend prior 
analytic calculations valid only in a few limiting regimes. Particular attention is given to the force 
direction, shown never to be directed opposite to the fiow except in the continuum limit, where 
other forces are of much stronger magnitude. 



Understanding the behaviour of a single dust parti- 
cle embedded in a flowing plasma is a physical prob- 
lem of practical relevance to the field of dusty plasmas, 
with concrete applications in fusion-grade, astrophys- 
ical, laboratory and technological plasmas [1]. The 
most important observable effect is perhaps momen- 
tum transfer between the plasma and the dust grain 
— the ion-drag force [2, 3, 4] — responsible for the dust 
dynamics as well as the formation of static configura- 
tions. Force calculations have very recently received re- 
newed attention in the context of weakly ionized plas- 
mas, where effects of fundamental physical importance 
transcending the field of dusty plasmas have been pro- 
posed. It has for example been proposed that two pos- 
itively charged dust particles can attract each other by 
a process analogous to cooper pairing in superconduc- 
tors [5]. 

More surprisingly, it is found that in the same 
strongly collisional limit the ion-drag on a single dust 
grain could be negative [6, 7, 8], implying that those 
particles see the surrounding plasma as a supcrfiuid [9] . 
Several published calculations have been done by treat- 
ing the plasma as a linear dielectric medium, and are 
therefore only valid for large shielding lengths As ^ 
690, where 690 is the 90° scattering parameter of the 
dust particle. Only the weak and moderate [3] as well 
as strongly collisional [6] regimes have been investi- 
gated with this approach, and the results of Ref. [3] 
are not directly quantifiable because they depend on 
the dust floating charge Qf, not known a priori. In 
the highly collisional limit, full non-linear calculations 
have been performed under the assumption of negligi- 
ble ion diffusivity (mobility regime) [7, 8], predicting 
that ion-drag reversal only occurs at strong coUisional- 
ity. Some Monte-carlo simulations [10] claim a negative 
ion-drag occurs also in the weak coUisionality regime. 

The purpose of this letter is to calculate the ion- 
drag force on a conducting spherical dust particle (ra- 
dius Tp) in the entire range of coUisionality, by solv- 
ing the full self-consistent non linear problem; this 
bridges the gaps between all available partial results. 
For this purpose we use the Particle In Cell (PIC) 
code SCEPTIC [4], recently upgraded to account 
for constant collision-frequency {u^) charge-exchange 
events [11]. More complex collision models could be 
implemented, but the present approach limits the num- 



ber of parameters while only negligibly affecting the 
physics [12]. SCEPTIC is a hybrid kinetic code where 
the coUisionless electrons have the Boltzmann density 
ne = TLoo exp(0) (0 = e^/Te and Uoo is the elec- 
tron density at infinity) and the ions are advanced ac- 
cording to Newton's equation between two collisions 
(rriidv/dt = —ZeS/V) up to convergence of the simu- 
lation. The electrostatic potential is solved on a two- 
dimensional spherical mesh centered on the dust parti- 
cle and extending several electron Debye lengths A^je ~ 
■^/eoTe/rioce^. The dust potential can be prespecified 
or floating, i.e. self-consistently set to a value such as 
to equate the ion current with the analytic coUisionless 
electron current /g = 4,Trr'ivte/2^/TT (vte = y^T^T/rTi^ is 
the electron thermal speed). At the outer boundary we 
impose dlncp/dr = -~rp{\/{Xsr) + 1/?"^), a condition 
valid at weak and moderate collision frequency as well 
as in the continuum limit if A^ 3> rp. 

The ions are reinjected at the outer boundary 
according to their distribution function at infinity 
("SCEPTICl" convention in Ref. [11]). We focus on 
the situation where the neutral background is station- 
ary and the ion flow driven by an external electric 
field Eext = "v^miVc/Ze (Z is the ion charge, usually 
Z=-|-l), whose role is to compensate the ion-neutral 
friction at infinity. In this case, provided vu ^ Vc {va 
is the ion- ion Coulomb collision frequency), the ion dis- 
tribution at infinity is given by [13]: 
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where erfcx(a;) ~ exp(a;^)erfc(x). f°° does not depend 
on i^c, and tends to a drifting Maxwellian /j\/ with 
temperature Ti and thermal speed vu = y^2Ti/mi as 
Vd/vu — *■ 0. Here Ti refers to the temperature of the 
neutral background; the effective ion temperature is 
T^^=T,+ m,v\ > T, and 7^^ = T,. 

The ion-drag Fi is the sum of the electrostatic force 
on the dust particle surface Fg arising from the in- 
teraction of its usually negative charge with the flow- 
induced anisotropy of the plasma (calculated by in- 
tegration of the Maxwell stress tensor at the particle 
surface) , the momentum collected by direct ion impact 
FP^, and the electron-pressure force FP = (averag- 
ing to zero on a conducting body [4]): F, = FP^ + Fg. 
Momentum conservation implies that in steady state. 



the ion-drag be also equal to the rate of momen- 
tum flux across any control surface surrounding the 
dust particle, in particular the outer boundary of the 
computational domain [4]. In a coUisionless plasma: 
Fi = F?^ + F^ + F° (F?^: net ion momentum flux 
into the computational domain; Fg and F°: integrals 
of the electrostatic stress and electron pressure on the 
boundary). Of course because the SCEPTIC domain is 
not infinite, the electrostatic stress and electron pres- 
sures at the outer boundary are not negligible. In the 
presence of ion-neutral collisions, one must also con- 
sider ion friction with the neutrals and the momentum 
provided by the external electric field inside the control 
volume: 



^ 150 

"g 100 






CO 
Q 
O 



-Total 
-E-field 
- Electrons 
Ions 
-Collisions 




10 H/o 10 

v^ / (ZT lv(\y'^ 
d ^ e r 

FIG. 1: Ion-drag force as a function of drift velocity com- 

7i,(x)(vd - v(x))d^x, (2) Puted by SCEPTIC for T. = 4.33 • lO^^ZTe, Xoe = 64.8rp, 

Domain (tip ~ —1.49 (Here imposed rather than floating), and 



where v(x) is the ion fluid (local average) velocity. 
F° can be either positive or negative; the integral in 
Eq. (2) is however convergent and F° tends to a limit 
as the domain size is increased. We will refer to Fim, 
Fe, Fe, F„ respectively as the "Ion", "E-field", "Elec- 
trons" and "Collisions" forces. 

We begin by comparing our code with the results of 
Schweigert and coauthors [10], who propose a weak 
collisionality regime {Xmfp = \/oTi/Trmi/iyc ^ A^) 
for which they find a negative ion-drag. The dimen- 
sional parameters corresponding to their Fig. (2) are 
Tp = A.l^im, As = 20 - lOO^im, Qp = 3.6 • lO'^e (dust 
charge), T^ = 6eV, T, = 0.026ey, P = 75 - 150Pa 
(neutral pressure), CTc = 3.53 • lO^^^cm^ (constant 
collision cross-section), rioo — 3.6 • lO^cm"'^. He- 
lium ions are used, but this information is irrele- 
vant because their work is based on a direct orbit- 
integration approach, where the dust charge as well 
as the potential distribution are prespecified. If we 
use the usual linearized formula for the sphere ca- 
pacitance: C = 4:TTeQrp{l + Tp/Xs), for simplicity 
set Aj^ = Xjj^ = A^^g -I- A^)^, and convert constant 
cross-section into constant collision-frequency accord- 
ing to Vc = o'c\/8Ti/TTmiP/Ti, we find the correspond- 
ing parameters to use for our SCEPTIC simulations. 
Fig. (1) shows the total ion-drag as well as its differ- 
ent components for the case corresponding to the curve 
As = 20fim and P = IbOPa in Fig. (2) of Schwcigert's 
pubhcation [10]. The ion-drag calculated by SCEPTIC 
at the dust surface and the outer boundary of the com- 
putational domain agree (their balance is of course dif- 
ferent), which gives us strong confidence that our code 
performs properly and the runs are well converged. 
More important the total ion-drag is positive, and does 
not agree with Schwcigert's according to which it can 
reverse for low drift velocities. 

Our hypothesis to explain Schwcigert's erroneous 
conclusion is as follows. Fig. (2) shows the two ion- 
drag curves from Fig. (2) of Ref. [10] with As = 
20 i^im (lines), as well as the ion-drag computed by 
SCEPTIC at the outer boundary of the domain with- 
out accounting for F° (lines with circle-markers), 
and the correct ion-drag (lines with square-markers). 



fc ~ 6.27 • 10 '^y''ZTe/mi/rp; corresponding to the curve 
Xs = 20/im and P = 150Pa in Fig. (2) from Ref. [10]. 
"Inner" (solid lines) and "Outer" (dotted lines) refer to 
the forces calculated at the dust surface and outer bound- 
ary of the computational domain; of course the total force 
(squares) does not depend on where it is evaluated. 



The strong similarity between Schwcigert's results and 
SCEPTIC-computed drags not accounting for F° sug- 
gests that the far-from-negligible F° has been ne- 
glected in Ref. [10]. 
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FIG. 2: Ion-drags from Fig. (2) in Ref. [10] at A^ = 
20/im ("Schw."), ion-drags computed by SCEPTIC at 
the outer boundary of the computational domain ("To- 
tal"), and ion-drags computed by SCEPTIC without ac- 
counting for F° ("Total- Collisions"). Pi: P = 75Pa 
i.e. i^c = 3.14 ■ IQ-'"^ y^ZTe/nii/rp; P2: P = 150Pa i.e. 
Vc = 6.27 • 1Q~^ y/ZTTfrni/rp (P2 curves correspond to the 
parameters of Fig. (1)). The correct ion-drag is is always 
positive (lines with square-markers). 

The only regime where we have observed a nega- 
tive ion-drag is when Vc 3> ^J ZTe/rrii/rp. In this 
highly coUisional limit a fluid treatment is appropriate, 
and analytic calculations assuming negligible diffusiv- 
ity (i.e. mobility-dominated physics) and As 3> rp have 
recently been published [7, 8] . According to this model 
the ion density is uniform on the ion fluid stream- 
lines, hence the ion current (required for an analytic 
calculation of the floating potential) and the electro- 



static stress at the particle surface only depend on the 
stream-lines' topology [7, 8]: 
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where oq = ao/j^p = max fl, ^JrpmiVcVd/{Ze\Vp\) 
F^ is always negative or null, because in the considered 
regime the downstream shadow of the dust particle is 
fully ion-depleted [6] . The ion impact part of the ion- 
drag, depending on both the stream-lines topology and 
the ion dynamics along those lines, can be decomposed 
in two parts. The first corresponds to the flux of z- 
momentum to the dust if the ions were indeed cold [7] : 






~al 



Tal 



(5) 

The second accounts for the ion pressure at the dust 
surface, due to the effective ion temperature along 
the stream-lines, never negligible even if Ti <^ T,,: 
Ti ||(x) = Ti + TOiv(x)2 (The heating depends on the 
local fluid velocity v(x) rather than vj, because the 
stream-lines are curved in the dust vicinity). The 
ion temperature perpendicular to the stream-lines is 
Tj^_L = Ti. In the limit Ti ^ T^, integration of the ion 
pressure at the dust surface yields: 
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where a = —ZeVplimiVdVcfp) (we recall that ZeVp < 
0); FL^ and FL^ are positive. Eqs. (3,4,5,6) corre- 
spond to the mobility model. 

Fig. (3) is a plot of forces and floating potential 
versus coUisionality for H'^ ions with Ti = 0. 01 ZT^, 
that we now analyze from high to low coUisionality. 
In the continuum limit the mobility model (dashed 
lines) agrees with SCEPTIC results (symbols), thus 
validating both works. As the coUisionality decreases, 
the Uoating potential rises to a maximum at Uc — 
O.Oly^T^/rrii/rp (a moderate coUisionality tends to in- 
crease the ion current [11]). The E- field part of the ion- 
drag follows a similar trend, i.e. peaks when Xmfp ^ As 
(shielding length); the reason being that collisions 
favour ion focusing downstream by breaking orbital 
angular momentum conservation, hence increase the 
polarization field in the drift direction. This effect, 
carefully explained in Ref. [3] , is not taken into account 
by Maiorov [14] who therefore arrives at the same er- 
roneous conclusion as Schweigert (i.e. that a negative 
ion-drag can occur at weak coUisionality). The dip 
in ion-impact force seen at I'c — O.Qiy/T^Jrrii/rp is 
due to the local Uoating potential maximum, whose 
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FIG. 3: Ion-drag and dust floating potential against coUi- 
sionality for H+ ions; T, = O.OlZTe, Vd = ^ZTe/nn and 
Xoe = 20rp. The points correspond to SCEPTIC data at 
the dust surface, and the dashed lines to the mobility cal- 
culations (4,5,6,3). 



effect is to reduce the energy at which ions are col- 
lected. The local minimum in the E-field part of the 
ion-drag at Vc '^ 10"^ yj ZTe/mi/rp corresponds to 
Xrnfp ~ As. As we keep reducing the coUisionality, 
the floating potential decreases to its "coUisionless" 
{vii -^ Vc "^ \/ ZTe/nii/rp) value, causing the E-ficld 
part of the ion-drag to increase. 

Let us discuss further our results in a regime where 
the conditions 690 ^ As and Vd <C vu are satisfied. 
This allows us to directly compare our computations 
with the result of Ivlev and coauthors, who calculate 
the electrostatic part of the ion-drag using the lin- 
earized plasma equations, not accounting for down- 
stream depletion due to the finite-sized grain (appro- 
priate at weak or moderate coUisionality) [3] : 
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7r/8Am/p, InA = ln(AD/2&go), where 

The function K is given by Kix) = 
^-l)a;V(l + a;2)-y^ln(H-a;2). 

We also used the vacuum expression for the dust par- 
ticle capacitance: Vp = Qp/Airearp. 

Fig. (4) is a plot of forces and floating potential ver- 
sus coUisionality for equithermal Ar^ ions and elec- 
trons, Xoe = 30rp and Vd = 0.2y^ZTe/mi. The ion- 
drag components and dust floating potential follow a 
trend similar to what has been observed in Fig. (3), 
although perhaps with a more monotonic behaviour 
because at high temperature the ion flow is less sensi- 
tive to slight variations in floating potential. Of course 
the ion-drag peak is located at a higher coUisionality 
{Xmfp (x \/Ti is larger here). The agreement between 
SCEPTIC and Eq. (7) is only qualitative; recall that 
Eq. (7) is only valid to logarithmic accuracy in A^i/fego- 

Fig. (4) also shows the high coUisionality limit of (^j 
and Fe- While the mobility model does not apply for 
Ti = ZTf^ (It requires Ti <^ ZTg), it is valid in the 
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FIG. 4: Ion-drag and floating potential versus coUisionality 
for Ar+ ions; Ti = ZTe, va = 0.2i/Z7Vm~ and Xoe = 
30rp. The points correspond to SCEPTIC data at the dust 
surface, the dotted line to Eq. (7), and the dashed lines to 
the hmits (8,9). Foxt is the background electric force on 
the dust particle (Eq. (11)), directed upstream. 



limit S-o 3> 1, i.e. when E^xt strongly dominates the 
dust Coulombic field. In other words as do —^ oo: 



We have also performed extensive simulations, not 
shown here, where the ion drift is driven by a neutral 
flow assumed unaffected by the dust particle (kinetic 
regime); the ion distribution at infinity is then given 
by a drifting Maxwellian rather than by Eq. (1). The 
ion-drag is very similar in magnitude to the results 
presented here, for comparable coUisionality and flow 
velocity. The negative i^ext (Eq. (11)) is however re- 
placed by a positive neutral drag. 

In conclusion, we have presented fully self-consistent 
calculations of the ion-drag force on a spherical dust 
particle over the entire range of charge-exchange col- 
lisionality, when the ion flow is driven by an external 
electric field. We have shown that the ion-drag behaves 
similarly to the floating potential; it shows a local max- 
imum when Xmfp ~ As, and tends to a limit at high 
i^c- Although there are physically meaningful parame- 
ters where the ion-drag is negative, that occurs only in 
collisional regimes where the direct electric field force 
far exceeds the drag. 

Leonardo Patacchini was supported in part by 
NSF/DOE Grant No. DE-FG02-06ER54891. 
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regardless of Ti/ZT^ (provided of course the formula 
for vacuum capacitance is applicable, i.e. Xjje 2> rp). 
This proves that past a certain coUisionality, increasing 
Pc does not increase the negativity of Fe, bounded by: 
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where the logarithm is the floating potential calculated 
by equating Eq. (8) with the coUisionlcss electron cur- 
rent. 

We have therefore shown that only in the continuum 
limit (i^c ^ \/ZTe/ini/rp) can the ion-drag reverse. 
There the ambipolar electric force on the dust particle 

(Fext = QpEext) 
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(11) 
is much larger than ion-drag itself (Fig. (4)); in the 

typical situation where Qp ~ ^irtorpVp < 0, Fext is 
directed upstream. 
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